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ABSTRACT 

We present the first entirely ground-based astrometric determination of the proper motion 
for the Fornax Local Group Dwarf Spheroidal satellite galaxy of the Milky Way, using CCD data 
acquired with the ESO 3.5 m NTT telescope at La Silla Observatory in Chile. Our unweighted 
mean from five Quasar fields in the background of Fornax, used as fiducial reference points, leads 
to fiaCOsS = 0.62 ± 0.16 mas y""'^, and fig — —0.53 ± 0.15 mas y~^. A detailed comparison 
with all previous measurements of this quantity seems to imply that there is still no convinc- 
ing convergence to a single value, perhaps indicating the existence of unnacounted systematic 
effects in (some of) these measurements. From all available proper motion and radial velocity 
measurements for Fornax, we compute Fornax's orbital parameters and their uncertainty using 
a realistic Galactic potential and a Monte Carlo simulation. Properties of the derived orbits are 
then compared to main star formation episodes in the history of Fornax. 

All published proper motion values imply that Fornax has recently (200-300 Myr ago) ap- 
proached perigalacticon at a distance of '^150 kpc. However, the derived period exhibits a large 
scatter, as does the apogalacticon. Our orbit, being the most energetic, implies a very large 
apogalactic distance of ^ 950 kpc. If this were the case, then Fornax would be a representative 
of an hypervelocity MW satellite in late infall. 

Subject headings: Astrometry: Proper motions — Galaxies: Local group — Galaxies: Dwarf spheroidal 
galaxies — Galaxies: Fornax — Galaxies: Proper rriotion — Stars: proper motions 



1. Introduction 

Studying the kinematics of the satehites of the 
Milky Way (MW) allows us to adddress various 
fundamental issues such as the origin and evolu- 
tion of this satellite system (and the MW itself), 
the role of tidal interactions in the evolution of 
the Local Group (LG) and the matter distribution 
of the latter (including that of the dark matter, 
thus allowing tests of so me cosmological predic- 
tions ( Shava et al.ll2009l )). This requires tracing 
their positions back in time, by integrating their 
orbits, which, in turn, requires knowing their cur- 
rent positions and their full 3-D space velocities. 
While radial velocities for LG galaxies are known 
to better than '^5 km s^^, and their distances to 
~10%, the biggest source of uncertainty rests on 
their proper motions (PMs). 

In the year 2000 we started a ground-based 
program aimed at determining the absolute PM 
of three southern dwarf Spheroidal (dSph) galax- 
ies. Carina, Fornax and Sculptor, with respect to 
background Quasars (QSOs), used as inertial ref- 
erence points. Three to four epochs of homoge- 
neous CCD data were obtained using a single tele- 
scope-(-detector setup over a period of eight years: 
Ours is the first entirely optical CCD/ground- 
based PM study of an external galaxy other than 
the Magellanic Clouds. 

In Mendez et al. (2010, hereafter Paper I), we 
presented a detailed description of our methods, 
as well as our first results for the PM of Fornax, 
based on one QSO field (QJ0240-3434B). In this 
paper we present our final PM for Fornax based 
on measurements from five QSO fields in the back- 
ground of Fornax. In Section [2] we present a sum- 
mary description of our observational material, in 
Section [3] we explain how we obtained our PMs, 
and, finally, in Section 13.11 we compare our results 
to previous studies and present our main conclu- 
sions. 

The Fornax dSph galaxy, at a distance of 
147 kpc (e.g., Pietrzyhski et al. 2009), is relatively 
isolated, luminous and well-resolved into individ- 
ual stars. It seems to be dark matter dominated, 
and has an estimated total mass of ^ 10^ Mq (e.g.. 
Walker et al. 2007). Unlike other dwarf galaxies, 
Fornax harbours five globular clusters ([Hodge I 



and t appears to have a complex stellar 
substructure in the form of shell-like features in- 
dicative of recent merger activity (Coleman et al. 
2004, Coleman & Da Costa 200 5, Olszewski et 
al. 20 06). Furthermore, results bv lBattaglia et al 



(I2OO6I ) suggest that the ancient stellar population 
in the centre of Fornax is not in equilibrium (ap- 
parent as a non-Gaussian, double peaked velocity 
distribution), which also points to a relatively re- 
cent accretion of external material, such as gas 
accretion due to the merger with another smaller 
stellar system. 

2. Observations 

All observations were carried out with the "Su- 
perb Seeing Imager", SuSI2, attached to one of 
the Nasmyth focii of the ESO 3.5 m NTT tele- 
scope at La Silla Observatory in Chile. The overall 
ch aracteristics of th e det ector are fully describe d 



D'Odorico I (Il998l) and lD'Odorico et al.l (Il998l) 



We followed exactly the same observational, re- 
duction, and calibration procedures for all of our 
QSO fields as described in Paper I. Therefore, in 
as much as possible, our dataset is homogeneous 
from this point of view (for certain limitations to 
this statement, see Section [3]). SuSI2 is a mosaic 
of two 2kx4k EEV 44-82 chips (called #45, on 
the West side, and # 46, on the East side, for the 
adopted rotator angle of 0°). As explained in Pa- 
per I, we always placed our QSOs at a nominal 
pixel position near {x,y) = (3 000, 2 100), close to 
the middle of chip #46, and only data from this 
chip was subsequently used for our astrometry. 

All of the astrometric observations, including 
those required to compute the differential chro- 
matic constants (DCR - see Paper I, Section 3.3.), 
were acquired through a Bessel R filter, whereas 
for the blue frames needed to construct color- 
magnitude diagrams (CMDs) we used the Bessel 
B filter, both of which are part of the standard set 
of filters available for SuSI2. 

Our initial list of QSOs in the background of 
Fornax comprised all eleven distinct QSOs (one, 
reported in Paper I, forms a gravitational lens 



^Another notable exception is the tidally disrupted Sagittar- 
ius dwarf which, along with Fornax, are the most massive of 
the MW satellite galaxies after the LMC and SMC. Sagit- 
tarius appears to have globular clusters that are, or were, 
associated with it (e.g., Law & Majewski 2010). 
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pair) repor ted by Tinney et al.l ( 19971 ) (their Ta- 
ble 5) and [Tinnesl (ll999|) (his Table 1). Unfor- 
tunately, several of them proyed to be useless for 
astrometry: deep, good-seeing (FWHM^O.5 arc- 
sec) images taken during the first epoch revealed 
that they were either too faint for our required as- 
trometric S/N (> 200 integrated over the PSF fit- 
ting radius; see Paper I, Section 3.2 and Figure 4), 
had a noticeably elongated structure (and there- 
fore had a different, usually more extended, PSF 
than that of the stars) , had a very nearby (usually 
stellar-like) bright companion, or were definitely 
blended with field stars. As explained at length in 
Paper I, any of these features render these targets 
unsuitable for high-accuracy relative astrometry. 
These problematic QSOs were dropped from the 
observing list in subsequent epochs, and we con- 
centrated on "clean" (as far as we could determine 
with a seeing of ~0.5 arcsec) QSOs. Our final list 
of five QSOs, from which we were able to deter- 
mine the PM of Fornax, is given in Table [1] where 
we also show the full list of known QSOs behind 
Fornax from the above-cited references. 

In Figure [T] we show the stellar configuration 
in the immediate surroundings of each of our five 
selected QSOs. All of these images were acquired 
on the first epoch, when we consistently had very 
good seeing (FWHM'^0.5 arcsec). The only po- 
tentially problematic case was that of QJ0239- 
3420 which has a faint companion to the NW, at a 
distance of ^ 1.4 arcsec. Fortunately, this nearby 
source is very dim (with a Star/QSO peak bright- 
ness ratio smaller than 0.1), and far enough (far- 
ther than 1.2 x FWHM, adopted as the PSF fit- 
ting radius, see Paper I, Section 3.2 and Figure 4) 
that, even on our worse frames (with a FWHM of 
^1.0 arcsec), it did not pose a problem for the 
astrometric solution, and so it was fully included 
in our analysis below. 

Table [5] contains a summary of the observa- 
tional material that was, in the end, used to com- 
pute our PMs. We note that more data was ac- 
quired for these QSO fields, but proved to be use- 
less on account of bad seeing, poor image qual- 
ity (deteriorating the astrometric solution), or bad 
sky transparency. As explained in Paper I, at the 
start of our program we adjusted the exposure 
time on the basis of seeing (typical values were 
between 300 s and 900 s), but later, it was de- 
cided to use a fixed integration time of 900 s in all 



cases for simplicity. 

3. Proper motions 

To determine the PMs we used the same proce- 
dure described at length in Paper I. A flow-chart 
summary of the full process is presented in Fig- 
ure O We refer the reader to Paper I for further 
details on the methodology and precise meaning of 
all steps. Even though in our PM solutions we only 
included frames acquired within ±1.5 hr. from the 
meridian (see Table [2]), all of our coordinates were 
(pre)-corrected for (continuous) atmospheric re- 
fraction and DCR as described in Paper I. We also 
excluded from the PM solution all frames with a 
FWHM > 1.0 arcsec, as they clearly deteriorated 
the linear fit of Baricentric position vs. epoch di- 
agram (see, e.g., Fig. 12 in Paper I), where these 
frames stand-out due to their large scatter. 

All PM data were treated as homogenously as 
possible, including the following constraints: 

• All stars with /i > 2.0 mas y"^ in our ini- 
tial local reference system (LRS, basically a 
set of "high-quality" reference stellar images, 
eventually bona-fide Fornax stars - see Pa- 
per I), indicating that they are either fore- 
ground Galactic stars, or that they have a 
high (pseudo)-PM, possibly due, e.g., to a 
faint unresolved companion or another prob- 
lem in the image), were purged from the 
LRS. The full geometric registration (and 
PMs) for the remaining LRS stars and the 
QSO were re-computed in an iterative pro- 
cess, 

• All LRS stars that exhibited a registration 
residual > 3(7 of the formal rms of the two- 
dimensional geometric registration in the X 
or Y coordinates in at least four (not neces- 
sarily consecutive) frames, were eliminated 
from the LRS. For registration we used a full 
grd Qj-^Qj. polynomial fit, which has been jus- 
tified in Paper I, 

• All LRS stars with cr^ > cr^^ mas y~^ were 
excluded from the LRS. The cut value, cr^^ , 
was estimated for each field based on the dis- 
tribution of PM errors vs. magnitude (see 
Table O, 



3 



• LRS stars exhibiting a PM error larger than 
that of the bulk of the LRS stars at a given 
magnitude (even if they, individually, had 
cT/j < CT/ic) were eliminated. This cut was 
done visually, on plots of PM error vs. mag- 
nitude for each QSO field, 

• Finally, we expunged objects that in the 
CMD appeared not to belong to Fornax. 
This photometric cleansing was done by first 
calibrating our instrumental photometry fol- 
lowing the procedure described in Paper I 
(Section 4), and then by comparing our re- 
sulting CMD with that of Stetson (1997), 
which defines the main features of the For- 
nax CMD. 

In Table [3] we show a summary of what results 
from applying the above criteria to the initial set 
of LRS stars for each of our five QSO fields. Af- 
ter applying all previously described cuts, we ver- 
ified that we ended with a uniform (X, Y) distri- 
bution of the LRS stars (see, e.g.. Figure 15 in 
Paper I), and with a reasonable distribution in 
magnitude and color. The resulting CMDs for 
the LRS stars and the respective QSO are shown 
in Figure [S] We note that in Paper I we cali- 
brated our photometry approximately by adopt- 
ing a color oi (B — R) ^ 1.3 for the red-clump 
from lStetsoiTI and the QSO blu e ma gnitude 

from t he works bv lTinnev et al.l (|l997l ) and lTinnev 
()l999l) . However, we noticed that, probably as re- 
sult of uncertainties in the blue photographic mag- 
nitudes for the QSOs (and/or possible QSO vari- 
ability), the ordinate in these figures had a large 
zero-point variation from field to field, and hence 
the red-clump did not fall at the same apparent 
R magnitude for the five fields, as is expected. 
We therefore decided to adopt a calibration based 
on the color of the red-clump (as before), but fix- 
ing, instead, the magnitude of the red-clump at 
R ^ 20.6, also from the photometry by Stetson 
(1997, his Figure 7). 

In Figures |4] and [5] we show, respectively, the 
barycentric position vs. epoch and the vector- 
point diagrams for our five fields, for the final LRS 
stars and the corresponding background QSO. 

After processing all of the QSO fields follow- 
ing the protocols described in the previous para- 
graphs, wc computed PMs for the five QSO fields 
presented in the upper part of Table [TJ The re- 



sults are summarized in Table |4l and are plotted 
in FigureH TableHgives the PM in RA and DEC 
as well as the overall rms of the linear fit of the 
barycentric coordinates vs. epoch. We recall that 
the slope of this linear fit gives us the PM of the 
Fornax stars with respect to the QSO, as explained 
in Paper I. Figure [6] clearly shows that our five 
measurements scatter more widely than would be 
expected from their error bars, and therefore com- 
puting a weighted mean does not seem appropri- 
ate for these data. Instead, and not knowing the 
source of the additional uncertainty affecting our 
proper motion measurements, we opted for com- 
puting an unweighted mean over the resulting val- 
ues for all five QSO fields. From Table S] we thus 
arrive at the following values for the PM of the 
Fornax dSph galaxy (and the standard deviation 
of the mean): fj,a cos 6 = 0.62 ± 0.16 mas y~^, and 
= —0.53 ± 0.15 mas y~^. 

3.1. Comparison to other studies 

Only two astrometric determinations of the PM 
for the Fornax dSph g a laxy a re available, namely 
that bv iDinescu et al. (2004), based on a combi- 
nation of ground-based plates and Hubble- WFPC 
data, and that based exclusively on HST data 
by Piatek et al. (2007, hereafter PI07), who 
present revised values to those reported earlier 
(p re-CCD Charge Tra nsfer Inefficieny corrections) 
in lPiatek et al.l ((2OO2I) . 

Figure [6] compares our results from individ- 
ual QSOs, to those obtained by PI07. We have 
three QSO fields in common with PI07; their 
values, along with our measurements, are given 
in Table [S] From this Table we see that the 
difference of the mean PMs derived from the 
QSOs in common between these two studies is 
less than Icr of our error in the mean in both, 
RA and DEC. On the other hand, the overall 
weighted mean value reported by PI07, based on 
4 QSOs (their Table 3) of ^i^cosS = 0.476 ± 
0.046 mas y-\ and ^5 = -0.360 ± 0.041 mas y^^ 
differs by about la from our overall unweighted- 
mean PM values. Even though our results are 
not statistically inconsistent with those of PI07, 
it is apparent from Figure [6] that our individ- 
ual results do exhibit a much larger scatter than 
their results (cr^^^ cosS = 0.36 mas y~^ and a^^ — 



0.34 mas y vs. 



0.06 mas y and 



<^fj.s = 0.09 mas y for PI07). We have no expla- 
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nation for these differences. 

The PM values pubUshed so far, along with our 
own value are Rive n in Table |6l More recently, 
Walker et al. have used a non-astrometric 

method called "perspective rotation" which is 
also included in the table (for details about this 
method see Paper I), because it provides a com- 
pletely independent measurement of the tranverse 
motion of Fornax (albeit with a larger error than 
the more recent astrometric determinations). All 
of these values are plotted, for comparison, in 
Figure [T] From the plot and figure we can see 
that, in general, there is a good agreement be- 
tween all measurements; indeed, none of them 
depart by more than 2a from the straight av- 
erage of the values in Table IH with our result 
being however the most extreme in this sense. 
Also, averaging the results from all of the au- 
thors, we find that the PM in DEC shows a larger 
scatter (cr^^ — 0.16 mas y^^) than that in RA 
(c/jccoss — 0.064 mas y~^). This would suggest 
that (some of) these measurements are affected 
by unaccounted systematic effects. Of the obvi- 
ous culprits, we can mention that ground-based 
astrometric data are affected by DCR. In the 
case of Fornax data (DEC^^ —34°) acquired from 
the Southern Hemisphere, DCR should however 
mostly affect RA PMs (we stress that our data 
have been corrected for this effect as far as possi- 
ble, see Paper I). 

On the other hand, HST data, while not af- 
fected by DCR, can be affected by "Charge Trans- 
fer Inefficiency" (CTI) in the CCD detectors due 
to the very low sky background which is insuf- 
ficient to fill-in the empty charge traps on the 
detector. These traps evolve in time because 
they are produced by in-fiight radiation damage, 
and induce systematic (time dependent) position 
shifts in the d e tected sources (for details see, e.g., 
Bristow et al. ( 20061) . especially his Figure 10). 
HST data has also been corrected for this effect, 
again, as best as they c ould (compare, e.g.. Ta- 
ble 3 from lPiatek et al.l (|2002l ) (CTI-uncorrected) 
and PI07 (CTI-corrected)). We note however that, 
since the HST cameras were not oriented exactly 
along DEC in the parallel readout direction, which 
is the direction mostly affected by CTI (see, e.g.. 
Figures 2, 3 and 4 in PI07), it is difficult to as- 
cribe the DEC scatter to this problem. Also, as 
mentioned previously, HST data shows a much 



smaller intrinsic scatter than our measurements 
(Figure [S]), thus suggesting small remaining sys- 
tematics due to CTI. 

4. Discussion and outlook 

One of the ultimate goals of these studies is a 
determination of a reliable orbit for Fornax. It 
is interesting to compare the impact of the dif- 
ferent values given in Table IHl on the kinematical 
and orbital parameters that can be derived from 
these measurements. Table [7] gives the Heliocen- 
tric PMs and the corresponding tangential veloci- 
ties in Galactic coordinates for the different values 
reported in Table HI From this table, it is clear 
that our PM value will render one of the most en- 
ergetic orbits yet derived for Fornax. 

For illustration purposes, and also as a key in- 
gredient to compute and interpret the orbit of For- 
nax from the above motions in Galactic coordi- 
nates, we have computed velocities in the Helio- 
centric (HC) and Galactocentric (GC) reference 
systems. For a detailed description of these vari- 
ous reference systems, and th e equations relatin g 
them, the reader is referred to lCosta et al. I (l2009l) . 
They are shown in Table [51 along with model- 
independent kinematical and orbital parameters 
derived from the PM values taken from different 

studies. Given the large values of Vt „ „ in compar- 

GC 

ison with the derived from all PM measure- 

GC 

ments, it is clear that Fornax must be close to peri- 
galaticon at this time. Also, given the rather small 
ratio of Lz/L it is clear that there must be signifi- 
cant excursions of Fornax away from the Galactic 
plane, regardless of the adopted PM values. We 
note again that, of all measurements available, our 
PM measurement yields the most extreme orbit. 

We have computed Galactic orbits for For- 
nax by integrating back in time the equation of 
motion under a realistic three-component (disk, 
halo and spheroid) model Galactic potential 
( Johnston et al.l[l995[ ). and using the current po- 
sition and velocity (and their uncertainties, in a 
Monte Carlo scheme) as initial conditions. De- 
tails of the integrator, and the Monte Ca r lo sim - 
ulations, are discussed in iDinescu et al.l (|2004[ ). 



A particularly important feature of the integra- 
tor used is its care to conserve energy and to- 
tal angular momentum, obvious features that are 
howevever sometimes tricky to achieve over the 
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full orbit, particularly when the number of in- 
tegration steps ('^several thousands in our case) 
is large. The adopted gravitational potential is 
strictly axi-symmetric, and it does not include the 
Galactic bar, nor the M W spiral pattern. R ecent 
numerical simulations by I Allen et al.l (|2008l ). that 
include these non axi-symmetric components, in- 
dicate however that their effect, in particular in 
the orbits of Galactic globular clusters (and there- 
fore also in the case of the more external satellite 
galaxies), is minimal, and it should not alter our 
conclusions importantly. The results of these in- 
tegrations are shown in Table [9] and in Figure \8\ 
In Table |9] the meaning of the different columns 
is rather obvious, except perhaps for the last two 
columns that correspond to the crossing (vertical) 
velocity of Fornax through the Galactic plane (Vz) 
and the total speed at perigalacticon (Vp). As it 
can be seen, regardless of what PM values one 
adopts, comparing the Vp values to the Vt„„ in 
Table [8] confirms that the current Fornax position 
lies indeed near perigalacticon (see also Figure [HI 
top panel). Actually, all PMs indicate that the 
minimum distance, projected into the Galactic 
plane, happened 200 to 300 Myr ago (Figure [HI 
bottom pan el). This is a very in teresting quan- 
tity, because lBattaglia et al. I (l2006l) have found ev- 
idence of at least three distinct stellar components 
in Fornax: a young population (few 100 Myr old) 
concentrated in the centre of the galaxy, visible as 
a Main Sequence in the CMD; an intermediate age 
population (2-8 Gyr old, possibly related to a shell 
structure in Fornax, described in the next para- 
graph); and a n ancient populat ion (> 10 Gyr). 
More recently, iKirbv et al.l (|2011 ) have also found 
(see, e.g., their Figure 2) evidence of enhanced 
star formation in the range 200-300 Myr. One 
could therefore conclude that the latest episode of 
star formation on Fornax may have been indeed 
triggered by its perigalacticon passage. 

When tracing the galaxy back in time, the un- 
certainty in the computed position increases as 
time goes on, for a given uncertainty in the initial 
conditions, and, as we extrapolate further back in 
time, different PM values lead to quite different 
orbits. This is clearly shown in Figure [HI which 
shows the distance of Fornax from the Galactic 
center (^qq, top panel) and the same distance 
projected onto the Galactic plane {Rp\, bottom 
panel) as a function of time from now. Initially 



all PM values produce a similar dQ(j, R^^ vs. 
time, but the solutions diverge afterwards. As 
mentioned before, our value renders the most ex- 
treme solution, basically indicating that in a Hub- 
ble time Fornax has not completed an orbit yet. 
This is at odds with all of the other solutions 
which, while being different among themselves, do 
indicate nevertheless several perigalacticon pas- 
sages in the last 10 Gyr. If, as argued before, 
perigalacticon has had an influence on enhanc- 
ing star for mation , the stellar population results 
bylBattaglia et al.l ( 2006 ) favour a rather long or- 
bital period, thus supporting our solution. We 
note here that the 2 Gyr population belonging to 
the shell structure found by, Coleman et al. (2004j). 
Coleman fc Da Costal (12005^ and [Olszewski et al 



2006 h has been interpreted by these authors 
the product of a merger with a smaller, gas-rich 
system, and may not bear any relation to the inter- 
action between Fornax and the MW. We are thus 
seemingly left with two significant star formation 
episodes, one at (200-300 Myr), close to perigalac- 
ticon, and another one at age > 10 Gyr, which 
again argues for a longer orbital period (interest- 
ingly, from Figure [HI we see that 2 Gyr mark 
the most recent ap ogalactic positio n for Fornax 
for the PI07 and Walker et al.l (|2008[ ) orbits). Our 
extended orbit implies that Fornax would belong 
to one of the "hypervelocity" satellite s of the MW,^ 
as arg u ed by Kallivayalil et al.l ( 20061 ) . iBesla et al 



(|2007t ). IPiatek et all (|2008[ ) results all based 



the same HST data set) . One could perhaps won- 
der that the late infall nature of Fornax's orbits 
(assuming our PMs) is possibly related to the fact 
that it is the only satellite dwarf galaxy of the MW 
(Mateo 1998), along with the spatially dissipated 
Sagittarius dwarf (Law & Majewski 2010), that 
harbour a globular cluster population. 

A final word of caution regarding the above dis- 
cussion is that, from Table IHl we see that cur- 
rent PM measurements for Fornax implies that 
the range of derived orbital parameters is quite 
broad, and it seems adventurous to extract strong 
conclusions from them. As already mentioned, of 
particular concern is the possible existence of yet 
unaccounted systematic effects in the PM mea- 
surements. Only high-accuracy future astrometric 
satellite measurements, like GAIA, with expected 
uncertainties of a few /larc-sec (average over many 
stars) could help resolve this issue. 
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Fig. 1. — Finding charts indicating the configu- 
ration and surroundings of our five QSO fields. 
The QSOs are indicated by an arrow. The field-of- 
vicw is ±250 pix around the QSO (approximately 
40 arcsec on a side), and the FWHM of all images 
is ~0.5 arcsec. With the exception of QJ0239- 
3420, all QSOs arc well isolated from surround- 
ing sources (see text). Our tests indicate that the 
small source to the NW of QJ0239-3420 does not 
affect the astrometry. To illustrate the effect of 
seeing, in the bottom right inset, we show a zoom 
of the area around QJ0239-3420 for best (top) and 
worse (bottom) FWHM frames for this field. 
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Fig. 2. — Flow diagram from raw stellar pixel co- 
ordinates to final PMs for the QSO and bona- 
fide Fornax reference stars, based on procedures 
detailed in Paper I. We note that there is an 
elaborate previous selection of "suitable" (good 
S/N, isolated) reference stars (potential members 
of Fornax) from which precise pixel coordinates 
are derived through PSF fitting. All these steps 
ocurr prior to the entry po int in this flo w diagram, 
and are fully explained in lCosta et al. and 
Paper I. 
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Fig. 3. — Instrumental CMDs for our five QSO 
fields. The QSO is idicated by an ellipse. For 
comparison purposes among the different fields, 
our PSF photometry has been approximately cal- 
ibrated by adopting a red-clump magnitude and 
color of i? ~ 20.6 and {B — R) ~ 1.3 respectively 
(from Stetson 1997). 
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Fig. 4. — QSO barycentric position with respect 
to the (final) LRS stars vs. epoch for our five QSO 
fields. The lines are (unweighted) fits to the data 
points, and the negative of the slope of these fits 
corresponds to the PM of Fornax, whose values 
are tabulated in Table ID 
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Fig. 5. — Vector-point diagrams for our five QSO 
fields, after having applied the selection criteria 
described in the text. The QSO is indicated by a 
solid dot with error bars {la uncertainty on the 
slope of the barycentric position vs. epoch fit). 
The (pseudo-) PM error of individual LRS stars is 
similar to that of the QSO, so the mean motion 
of the bulk of the LRS stars has a very small un- 
certainty, and our final error is dominated by the 
PM uncertainty of the QSO itself (see Costa et al. 
2011). 
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Fig. 6. — Individual Fornax PMs determined from 
our five QSO fields (from Table IH open squares). 
Filled squares are the HST PC2-I-STIS results 
from four QSOs by Piatek et al. 2007, their Ta- 
ble 3. Three of these QSOs are common to ours 
(see our Table [5]). 
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Fig. 7. — Mean Fornax PM and their la er- 
rors as determined by vario us authors: Filled 
dot from iDinescu et al.l (|2004l) . fille d square from 
Piatek et al.l (|2007t ). filled star from 'Walke r et all 
POOSI ). and open square this work (see Table [6]). 
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Fig. 8. — X-Y and X-Z projections of the or- 
bit for Fornax from PM values by differe nt au- 



( 20041) ■ dot- 



Walker et al 



thors: Dashed line i s [Dinescu et al. 
dashed IPiatek et all ([20071), full line 
(|2008l) . and dotted hue this work. The filled dot at 
iX,Y,Z) = (-41.4,-50.9,-133.9) kpc indicates 
the current position of Fornax. 
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Fig. 9. — Fornax distance (log scale) as a func- 
tion of age from our orbit integrations, for PM 
values from different authors. Top panel is for 
the Galactocentric distance, bottom panel is for 
the distance projected on the plan e of the MW 



ble dot-dashed line for 



Both panels: Dot-dash ed line is for [ Dinescu et al 
(|2004l ). dashed line for IPiatek et all (l2007l) . dou 



Walker et all (|2008D . and 



dotted line for this work. 
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Table 1 

QSOS IN THE BACKGROUND OF FORNAX, FROM ITiNNEY ET AlI (jl997l ) AND ITINNEyI (jl999h . IN ORDER 

OF INCREASING RA*^ . 



QSO ID 


RA (J2000) DEC 


B 




9'= 


Comments 




[hh;mm:ss.s] 


[":':"] 


[mag] 


[arcmin] 


[dog] 




QJ0238-3443 


02:38:43.8 


-34:43:53 


20.2 


22.9 


222 




QJ0238-3440 


02:38:55.6 


-34:40:45 


20.1 


19.0 


223 




QJ0239-3420 


02:39:49.0 


-34:20:00 


20.6 


07.3 


344 




QJ0240-3434B 


02:40:08.2'^ 


-34:34:22'= 


19.9 


07.6 


166 


Roportcd in Papor I 


QJ0240-3438 


02:40:38.7 


-34:38:58 


20.2 


14.5 


146 




QJ0239-3425 


02:39:32.9 


-34:25:25 


20.4 






No finding chart available 


QJ0239-3421 


02:39:36.9 


-34:21:30 


19.9 






Fuzzy, elongated QSO image, with a FWHM 48% larger than that of stars in the FOV 


QJ0240-3434A 


02:40:07.7'' 


-34:34:20'' 


19.1 






Blended to foreground star 


QJ0240-3437 


02:40:19.0 


-34:37:20 


17.9 






QSO exhibits complex structure 


QJ0241-3420 


02:41:57.9 


-34:20:49 


20.4 






No finding chart available, far from Fornax's center 


QJ0242-3426 


02:42:06.5 


-34:26:12 


21.8 






Faint and far from Fornax's center 


QJ0242-3424 


02:42:19.9 


-34:24:20 


21.5 






Faint and far from Fornax's center 



^Thc upper part of the table shows the QSOs used in this work. The bottom part indicates those which could not be used for astrometric purposes or were otherwise 
discarded for various reasons, as indicated under column "Comments". 

^To compute the angular d istance (p) and position angle (0) of the QSO with respect to Fornax, we have taken the J2000 coordinates for the center of Fornax 
from Table 1 of iMateo I (1199811 . namely RA^02:39:59, and DEC^-34:27.0. We note that these values differ by more than 1 arcmin in RA and 3 arcmin in DEC from 
-<I the values presented in the review of LG galaxies by van den Bergh (1999, Table 1) which, given the low surface brightness (~23.4 mag arcsec"'^ in the V band) and 

large angular extent (a core radius of 14 arcmin) of the galaxy (Mateo 1998), is probably understandable. 

According to lTinnev et ahl l]1997|) component B is 5.9 arcscc East and 1.6 arcsec S of component A, whose coordinates are given in ^. In our images we measure 
similar offset values for component B, namely, 5.8 arcscc East and 1.4 arcscc S from component A. 



Table 2 

ASTROMETRIC R-BAND^ FRAMES effectively USED IN EACH OF OUR FIVE QSO FIELDS. 



QSO ID 


Epoch range 


# epochs 


# astrometric frames'^ 


# DCR frames 


DCR HA range 




[year] 








[hour] 


Q,I0238-3443 


2000.98 ... 2008.82 


4 


09 


09 


0.68 ... 2.99 


QJ0238-3440 


2000.98 ... 2008.82 


4 


19 


13 


0.36 ... 3.86 


QJ0239-3420 


2000.99 ... 2007.87 


3 


14 


12 


0.55 ... 3.70 


QJ0240-3434B'^ 


2000.61 ... 2007.85 


3 


15 


13 


-0.65 ... -4.03 


QJ0240-3438 


2000.61 ... 2008.83 


4 


20 


10 


-0.74 ... -3.62 



^As explained in the text, B-band frames for building up CMDs of each field were also acquired. 

^All of these frames satisfy that \HA\ < 1.5 hr, and FWHM < 1.0 arcsec (as an example, see Table 1 of Paper I). 

^Reported in Paper I. 



Table 3 

Local reference stars and cuts applied in each QSO field. 



QSO ID 


# initial LRS stars 


# final LRS stars 


[mas y ^] 


PM error vs. mag 
# stars purged 


CMD cleansing 
# stars purged 


QJ0238-3443 


295 


226 


0.78 


3 


12 


QJ0238-3440 


217 


175 


0.60 


6 


7 


QJ0239-3420 


337 


250 


0.60 


12 


11 


QJ0240-3434B 


260 


217 


0.50 


11 


7 


QJ0240-3438 


156 


123 


0.60 
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Table 4: Fornax PM from our five QSO fields. 



QSO ID 


COS 5 


^J'5 


a fit in RA 


a fit in DEC 




[mas y"-^] 


[mas y"-^] 


[mas] 


[mas] 


QJ0238-3443 


0.90 ±0.15 


-0.80 ±0.12 


1.5 


1.2 


QJ0238-3440 


0.55 ±0.13 


-0.66 ±0.15 


1.6 


2.0 


QJ0239-3420 


0.94 ±0.17 


-0.78 ±0.27 


1.8 


2.8 


QJ0240-3434B 


0.64 ±0.08 


-0.01 ±0.11 


0.9 


1.3 


QJ0240-3438 


0.06 ±0.21 


-0.40 ± 0.20 


2.8 


2.7 


Uiiwciglit (h1 mean 


0.02 ± 0.16 


-U.33 ± 0.15 
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Table 5: Comparis on of Fornax PMs derived by 
Piatek et alj pOOTI ) with our results for the same 
QSO fields. 



QSO ID 


Ha cos S 
[mas y-^] 


HS 

[mas y""'^] 


fla COS S 

[mas y~-^] 


HS 

[mas y~-^] 




1 

Piatek et al 


1 

. f2007) 


1 

This work 


QJ0238-3443 
QJ0240-3434B 
QJ0240-3438 


0.45 ±0.07 - 
0.42 ±0.10 - 
0.54 ±0.16 - 


0.39 ±0.06 
0.48 ±0.11 
0.33 ±0.16 


0.90 ±0.15 
0.64 ±0.08 
0.06 ±0.21 


-0.80 ±0.12 
-0.01 ±0.11 
-0.40 ±0.20 


Unweighted mean 


0.47 


-0.40 


0.53 


-0.40 
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Table 6 

Published mean Fornax PMs and their uncertainties. 



Reference 


/^a cos 5 
[mas y"^] 


Ms 

[mas y"^] 


PM 

[mas y'"^] 


PA" 
[deg] 


Comments 


Dinescu et al. ('20041 
Piatek et al. (20071 
This work 


0.59 ± 0.16 
0.476 ± 0.046 
0.62 ± 0.16 


-0.15 ± 0.16 
-0.360 ± 0.041 
-0.53 ± 0.15 


0.61 ± 0.16 
0.597 ± 0.044 
0.82 ± 0.16 


104 ± 15 
127 ± 4 
125 ± 11 


Plates, 48 galaxies and 8 QSOs 
HST+PC2+STIS, 4 QSO fields 
NTT+SuSI2, 5 QSO fields 


Walker et al. (2008') 


0.48 ± 0.15 


-0.25 ± 0.14 


0.54 ± 0.15 


118 ± 15 


Radial velocities, "perspective rotation" 



'^Position angle measured from North to East. 



Table 7: Proper motions in the Galactic coor- 
dinate system^ and tangential velocities'^ derived 
from the PMs given in Table El 



Reference 






Vi 


Vb 


Vt 




[mas y^"'] 


[mas y^"""] 


[km s""'^] 


[km s-i] 


[km s-i] 


Dinescu et al. (2004) 


0.04 ±0.16 


0.60 ±0.16 


31 ± 112 


423 ± 112 


424 ± 112 


Piatek et al. C2007) 


0.271 ±0.041 


0.532 ±0.046 


189 ±31 


370 ± 34 


416 ±33 


This work 


0.41 ±0.15 


0.70 ±0.16 


288 ± 107 


490 ± 113 


568 ± 112 


Walker et al. (2008) 


0.16 ±0.14 


0.52 ±0.15 


113 ±98 


360 ± 105 


377 ± 104 



"To derive Galactic motions from the J2000 fia cos (5 and fis PMs we have ado pted a position for the J2000 Galactic Pole of 
RA=(12:51:26.2754, DEC= +27:07:41.705) following [Miyamoto fc Somal l|l993l 1 (see, especially, their equation (29)). 

''To derive the tangential velocities Vi and Vt, we have a adopted a distance modulus for Fornax of (m — M)o = 20.84±0.15 from 
the "Araucaria distance scale project" (Pietrzynski et al. 2009). The velocity values above include a full statistical propagation 
of PM and distance (taken to be d = 147 it 10 kpc) errors. Vt is the total Heliocentric tangential velocity, and its error. 
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Table 8 

Fornax current Heliocentric^ and Galactocentric*^ velocities and specific orbital 
kinematical parameters from values of different authors. 



R.cfcrcnce 


fHC 

[km s-^] 


[km s-'-] 


^HC 
[km s-i] 


[km s"i] 


^'GC, 
[km s-^] 


^^GC, 

[km s-^] 


^'GC, 
[km s^i] 


[xlO^ (km s"^)^] 


[xio" kpc km s"^ ] 


[xlO^ kpc km s"^ ] 


^incscu ct al. (2004) 
"j^iatck ct al. (2007) 
This work 


-195 ± 109 
-36 ± 31 
-13 ± 106 


-359 ± 105 
-404 ± 31 
-550 ± 104 


126 ± 46 
104 ± 15 
154 ± 47 


-23 ± 62 
-32 ± 19 
-33 ± 62 


59 ± 107 
-93 ± 32 
-203 ± 105 


257 ± 98 
189 ± 29 
298 ± 96 


263 ± 98 
210 ± 30 
361 ± 99 


35 
23 
66 


38 
31 
53 


-3.9 
6.1 
13.3 


Walker et al. (2008") 


-95 ± 97 


-355 ± 96 


100 ± 43 


-28 ± 58 


-16 ± 97 


186 ± 89 


186 ± 89 


18 


27 


1.1 



^Heliocentric velocities have been computed using the adopted J2000 RA and DEC coordinates for the Fornax center (see TablefTll, which correspond to Galactic coordinates {l-,b) - 
(237.10°, —65.65°). The errors correspond to the statitical propagation of all observational errors. The (U, V, W) is a right-handed system that points to the Galactic center, Galacti 
rotation and North Galactic Pole respectively. 

^Galactocentric velocities have been computed adopting a solar peculiar veloicity of (■uq , vq , wq ) — (10. 00 ±0.36, 5. 25 ±0.62, 7. 17ib 0.38) km s~ ^ from fPehnen & Binneyl jl998D , a 
LSR speed of 220 ± 5 km s~^ and _R0 ^ 8.5 ± 0.5 kpc. The errors correspond to the statitical propagation of all observational errors. While these are commonly used values, we not 
that recent studies seem to suggest a larger V-component of the solar motion with respect to the LSR, ~13 km s~^ (e.g., Co§kunoglu et al. 2011). However, giv en the current 
uncertainties, this change will have little effect on the calculated motion of Fornax. Also, note that our uncertainty for the LSR speed is rather optimistic, for example JBovy et al.| (12009 
quote Vlsr, ^ 244 ± 13 km s"^, while [Brunthaler et "all (120111) argue for l^LgR ^ '^^^ ± 7 km s~^. 

'^Kinetic energy per unit mass, as seen from a point at rest with respect to the Galactic Center. 



'^Total L-^ + Ly + L^) and z-component of the angular momentum per unit mass, with respect to the Galactic Center, in our right-handed system. 



Table 9 

Fornax orbital parameters from PM values of different authors'^ 



Reference 


Orbital period 


Apogalactic distance 


Pcrigalactic distance 


-^max 


7 

mm 


Exccntricity 


Inclination 


Vz 






[Gyr] 


[kpc] 


[kpc] 


[kpc] 


[kpc] 




[dcg] 


[km s"^] 


[km s"^] 


Dinescu et al. (2004) 


8 ± 4 


328 ± 226 


148 ± 10 


223 ± 84 


-238 ± 108 


0.38 ± 0.21 


68 ± 8 


189 ± 9 


265 ± 51 


Piatek et al. ('2007') 


6 ± 1 


197 ± 58 


142 ± 14 


163 ± 35 


-166 ± 37 


0.16 ± 0.09 


76 ± 4 


186 ± 2 


221 ± 17 


This work 


21 ± 8 


956 ± 376 


148 ± 8 


498 ± 191 


-411 ± 172 


0.73 ± 0.18 


41 ± 8 


181 ± 19 


363 ± 51 


Walker et al. (20081 


5 ± 3 


164 ± 109 


131 ± 30 


147 ± 73 


-148 ± 66 


0.11 ± 0.18 


83 ± 7 


191 ± 8 


212 ± 41 



^Thc errors were computed as the inner 50% interquartile range derived from 1,000 simulations with uncertainties drawn from Gaussian distributions of the errors in PMs, W 
and distance. The integration period was chosen to be large enough (50 Gyr) so that the longest orbital period solution completed at least two full orbital cycles. 



